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Abstract 

Within the Biittiker dephasing model, the backscattering in the dephasing 
process is ehminated by setting a proper boundary condition. Exphcit expres- 
sion is carried out for the effective total tunneling probability in the presence 
of multiple pure dephasing scatterers with partial coherence. The derived 
formula is illustrated analytically by various limiting cases, and numerically 
for its application in tunneling through multi-barrier systems. 
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To simulate the phase-breaking effect in partially coherent transport through a meso- 
scopic system, Biittiker proposed a conceptually simple model by coupling electronic reser- 
voirs to the conductor The dephasing reservoir can be thought of as either a fictitious or a 
real branch voltage probe. Although this approach appears to be purely phenomenological, 
it however can be justified from a microscopic theory with proper approximations,i^ by 
viewing that both the electron-phonon interactions and the dephasing reservoir can be de- 
scribed by a self-energy function. Owing to the simplicity, the Biittiker dephasing model has 
received noticeable attention,li@ and been applied to transport through various mesoscopic 
systems.liHii 

Noticeably, in the original work of Biittikeifli and the later applications mentioned above, 
in addition to randomizing the electronic phase, the phase-breaking scatterer would also ran- 
domize the electronic momentum. Randomization of momentum means backscattering in the 
dephasing process, thus introduces an additional resistance. This undesired feature has been 
noticed and analyzed by a few authors,i2F§l commonly following the idea by coupling two 
voltage probes to model a single pure dephasing scatterer. In this paper, based on the origi- 
nal work of Biittiker (i.e. using a single reservoir to model a single dephasing scatterer) ,0'i 
an explicit expression will be derived for the effective total tunneling probability through 
a multi-barrier mesoscopic system in the presence of multiple pure dephasing scatterers 
with arbitrary dephasing strength. The underlying physics and practical application will be 
illustrated clearly. 

In general, consider the tunneling through a series of barriers shown in Fig. 1, where the 
squares stand for tunnel barriers, and the triangles for dephasing scatterers. They can be 
described in terms of scattering matrices as follows. For the individual (symmetric) barrier 
(e.g. the jth one), the tunneling property is characterized byi 



(1) 



which connects the two incoming amplitudes with the two outgoing amplitudes as is well 
known, tj and are respectively the transmission and reflection coefficients through the 



barrier, and can be characterized by a real parameter (0 < 5j < 1) as tj = ySj, and 
rj = i^l — 6j. For each dephasing scatterer, following Biittiker's approachEli, a dephasing 
reservoir is coupled to the system via two channels. During tunneling process, the electron 
has certain quantum probability of being scattered into the reservoir, undergoing phase 
randomization in it, then returning into the system via the coupler. As a result of dephasing, 
the re-emitted component does not interfere with that having not entered the reservoir. As 
a specific model, the coupler is described by the following scattering matrixS'i 
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This 4x4 matrix connects the 4 incoming amplitudes with the 4 outgoing amplitudes along 
the attached 4 channels. Here, e G [0, 1] is a parameter characterizing the dephasing strength 
ranging from complete coherence to complete incoherence. 

To carry out the effective tunneling probability from channel 1 to channel 2 shown in 
Fig. 1, the (2A^ + 2) x (2A^ + 2) S-matrix S^^^ of the entire system is needed. Note that 
the dimension of S*-^-* corresponds to the total number of external channels. By applying 
either a direct matrix algebra or the Feynman path technique, S*-^-* can be carried out via 
the following recursive expressions:i0 
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Here J, J = 1, 3, 4, ■ ■ ■ , 2A^ - 1, 2A^, and i,j = 2, 2iV + 1, 2A^ + 2. S^^^^^ is the scattering 
matrix for the subsystem of tunnel barriers (correspondingly, A^ — 1 dephasing scatterers); 
and the 4x4 scattering matrix s^^-* is for the A^th single scatterer-barrier segment. Z^r = 



XSiT] \ resulting physically from the multiple reflections.! x denotes the plane- 



wave propagation-induced phase-change factor for the electron moving between the two 
composite subsystems. For simplicity, we hereafter set x = I5 which is valid for far off- 
resonance tunneling. 

In practice of computing the S-matrix of the entire system, we start with S^^-* for the 
simple system of barrier-scatterer-barrier. Applying the recursive technique described above 
and making use of Eqs. (|Tp and (^, we obtain 



ri + a r2 atit2 
atit2 r2 + a^ri 
Pti al3rit2 



aPtir2 



Pt2 



a(3rit2 
-a(l - rir2) 



aPtir2 

Pt2 

-a{l - rir2) 

(3^2 



(4) 



Here a = y/l — e and P = ^/e are introduced to simplify the notation, and Zi = (1 — 
a^rir2)^^. The most important advantage of the present technique is its ability in treating 
multiple dephasing scatterers. To this end, we need the S-matrix of the added segment 
of scatterer-barrier. It can be easily obtained from Eq. by setting ri = and ti = 
1, and replacing t2 (^2) by tj (rj) for the jth coherent barrier. We denote the resulting 
4x4 S-matrix by s'--'^ [c.f. Eq. (^)]. Using Eq. (|) and the recursive rules of Eq. (^, we 
can therefore evaluate the S-matrix S^^-* for the entire system, thus calculate the effective 
tunneling probability through the system (i.e. from channels 1 to 2). In the following, for 
the convenience of comparison, we first present result under the momentum randomization 
boundary condition. Then, the formal result and its implication/application based on the 
pure dephasing boundary condition are carried out comparatively. 

Momentum Randomization Boundary Condition. Under this boundary condition, the 
re-emitting current from the nth reservoir into the system is assumed to be equally injected 
through its {2n + l)th and {2n + 2)th channels. We denote this current as J„. In the 
dephasing process, the electron number is conserved, thus the net current of the two channels 
connecting the system and the reservoir is zero. This feature is characterized by 
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(2I-t)J5 = K^ Vl, (5) 

where I is the N x N unit matrix, and = (Ji, J2, ■ ■ ■ , Jn)'^- Hereafter the superscript 
(■ ■ means transposition of matrix. Jl is the current injected from the left channel 1. T 

is a X matrix with elements fnm = T2n+l,2m+l + T2n+l,2m+2 + T2n+2,2ra+l + 72n+2,2m+2, 

where Tjj = Tij = ISij]"^. K*'^'* is a x 1 matrix with elements K^^ = T2n+i,i + 72n+2,i- 
Substituting the formal solution of Eq. (|^), = (21 — T)^^K^^'' J^, into the expression of 

~ (2) ~ 

the outgoing current in the right channel 2, = T21 Jl + K J^, we obtain the effective 
tunneling probability under the momentum randomization boundary condition 

= Jr/Jl = T21 + K^'\2l - T)-'K^'^ , (6) 

~ (2) 

where K is an 1 x matrix with elements K^^ = T2^2n+i + ^2,2n+2- 

Backscattering-Free Boundary Condition. To elucidate this boundary condition, let 
us consider in more detail the scattering on the dephasing scatterer (e.g. the nth one), 
see Fig. 2. The scattering matrix of Eq. (Q) relates the outgoing amplitudes = 
(«2n+i>«2n+2)^2n+i'^2n+2)^ to the iucomlug amplitudes a„ = (a2„+i, a2n+2, &2n+i> ^2)^ in 
terms of = S^a„. It is easy to check that the unitary property of and its specific 
structure would lead to the conserving relations 

|Q'2n+lP + I^2n+2P = k2n+2 P + |&2n+lP , (7a) 
|Q'2n+2r + 1^271+1^ = |«2n+l|^ + |&2n+2|^ • (7b) 

Now we make use of the pure dephasing condition: |a2n+iP = k2n+2p? ^^d 1 02^+2 P = 
l'^2n+iP- physical meaning is quite clear by noting that the phase-breaking only random- 
izes the phases of the forward and backward going waves, but does not cause any reflections 
. Because of this nature, we term it as a backs cattering- free boundary condition, which 
straightforwardly leads to 
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With this physical insight, we are ready to derive the effective tunnehng probabihty under 
the present boundary condition. We denote J„ = {j2n+i, j2n+2)'^ for the injecting currents 
from the nth reservoir into the system via the (2n + l)th and (2n + 2)th channels, and 
'^'n = (J2n+i5 J2n+2)^ ^^e curreuts being scattered into the reservoir from the system via 
the same channels. The pure dephasing boundary condition, now denoted as j2n+i — J2n+2, 
and j2n+2 ~ J2n+i, cau be expressed compactly for all the N reservoirs as 

J'b = BJb, (9) 

where 3b = {Ji, -^2, ■ ' ' > Jn)'^, J'b = {J[, J2i' ' ' 1 J'n)^ 1 ^'^d 6 = 1® a. aisa2x2 matrix 
with elements aij = 1 — 5ij. Again, a simple current counting leads to: 

Jii = T2iJi + K(2)jB, (10a) 
J'b = K«Jl + TJb. (10b) 

Here, K^^'' is a 2N x 1 matrix with elements T„i, and K*-^^ an 1 x 2N matrix with elements 
T2ni with n = 3, 4, ■ ■ ■ , 2N + 2. T is a 2N x 2N matrix with elements T^n, with m,n = 
3, 4, ■ ■ • , 2A^+2. Straightforwardly, by substituting Eq. (^) into Eq. (|TU]) an elegant expression 
is obtained for the effective tunneling probability 

Teff = Jr/Jl = T21 + K(2)(B - T)-iK(i) . (11) 

This equation is formally similar to Eq. (^. Both Eqs. (^) and (|11]) contain a common 
coherent term (the first one), and an incoherent term (the second one). However, they 
physically differentiate from each other. Below we detail our discussions on their difference 
by focusing on the special case of one dephasing scatterer. 
In the case of iV = 1, from Eq. (P) we easily obtain 

Teff = T21 + (T23 + T24)(2 - Tll)~^(Ti3 + Tm) 

= T21 + S,Sf/{S, + Sf) . (12) 

Here, following Biittiker,0i the incoherent backward and forward scattering probabilities 
are introduced: Si, = T13 + T14, and Sj = T23 + T24. In deriving Eq. (0), the symmetry 
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Tmn = Tnm and the sum rule J2mTmn = 1 have been used. On the other hand, for = 1, 
Eq. (|Tl]) leads to 

TeS = T21 + [(T44T23T31 + T33T24T41) 

+ (1 - T34)(T23T4i + T24T3i)]/Z , (13) 

where Z = {1 - T^^y - T^sT^- 

Further, in the completely incoherent regime (e = 1), T13 = Ti = T24 = T2 = |t2p, 
^33 = 1 — ^17 and T44 = 1 — T2. Other Tmn in Eqs. ( p^ and ([131) are zero. In this regime, 
Eq. (112]) reduces to 



^eff 

whereas Eq. (|13D gives rise to 



1 1 



(14) 



eff 



n -1 



Ti +T2 



(15) 



The latter equation (0) can also be obtained by a rather simple treatment based on incoher- 
ent multiple reflections and transmissions through two barriers. The interesting difference 
between Eqs. (|T^) and (|T3p is highlighted as follows. It is well known that the two-terminal 
resistance is related to the effective total transmission probability Ttot (i-e. Tes and Tes) via 
the Laudauer formula,!! 7?,tot = (^/e^)(l/Ttot)- From Eq. (0) the system resistance can be 
expressed as 

ntot = 7^l + 7^2 + 7^, + 7^, . (le) 

In this decomposed form, TZj = {h/e^){\ — Tj)/Tj (j = 1,2) is the Landauer resistance 
for the jth conductor, TZc = h/e^ is the so-called contact resistance rooted in the two- 
terminal configuration (measurement ).iS Interestingly, the dephasing scatterer contributes 
a constant resistance TZg = h/e^ in the completely incoherent regime, due to the momentum 
randomization. On the other hand, from Eq. (|15D we have 



Titot = 7^l + 7^2 + 7^, . (17) 



We see that under the boundary condition Eq. (|^) the dephasing scatterer only plays a 
phase-breaking role. As a result, the Landauer resistances are connected in series in a 
purely classical way, and the dephasing source does not cause additional resistance. 

Interestingly, the backscattering-free nature on the dephasing scatterer can be further 
elucidated by considering the simple transmission through only a dephasing scatterer (i.e. 
with no tunnel barriers), with arbitrary dephasing strength. In this case, the relevant 
transmission coefficients are: T21 = T34 = 1 — e, T13 = T24 = e, T33 = T44 = T14 = 
T23 = 0. Accordingly, Eq. (|l^) gives rise to Teg = 1 — e/2, whereas Eq. ([T3|) leads to 
Teff = 1. This result clearly shows the distinct nature of the two dephasing models: one 
causes backscattering, another is backscattering free. 

For two resistors connected in partial coherence (0 < e < 1), the simple sum rule of the 
individual resistances as Eqs. (pIB]) and ([T7|) breaks down. However, the resistance difference 
ATZ = (/i/e^)(l/Tefj — 1/Teff) based on Eqs. (|1^) and (|13]) is a proper quantity to characterize 
the additional resistance caused by the dephasing scatterer under the the momentum ran- 
domization boundary condition. Figure 3 shows ATZ as a function of the dephasing strength 
e. In general, the interplay between the backscattering on the dephasing scatterer and the 
tunneling through the individual barriers leads to ATZ{e) depending on the barrier-tunneling 
strength as shown in Fig. 3. 



To further illustrate the application of Eqs. (H) and ([TTD in combination with the recursive 
rule of Eq. (0), we briefly present results for tunneling in the presence of multiple phase- 
breaking scatterers. For clarity. Figure 4 shows the relative tunneling probabilities versus 
the dephasing strength, i.e., Teflf(e)/Tcfj(e = 0) and Tefj(e)/Teff(e = 0), by sohd and dashed 
curves, respectively. In the weak tunneling regime {S = 0.1) shown in Fig. 4(a), dephasing 
enhances the tunneling remarkably, and the two dephasing models, i.e., Eqs. (^ and (^TJ), 
give almost the same results. However, in the strong tunneling regime {6 = 0.9) shown in 
Fig. 4(b), the two dephasing models give noticeably different results. More interestingly, the 
turnover behavior predicted by Eq. (^, which was paid special attention by Buttiker,0 does 
not occur in the pure dephasing model of Eq. (|ll|). Therefore, the turnover behavior appears 



a consequence of competition between the incoherent-tunneling and the dephasing-induced 
backscattering. It disappears after the backscattering being ehminated. 

In summary, we have presented a unified treatment for phase-breaking tunneling based 
on the Bttiker model. Simple adopting of two types of boundary conditions can clear up the 
distinct natures of two dephasing models. The derived explicit expressions in combination 
with the recursive rules can be conveniently applied to tunneling through multi-barrier 
systems in arbitrary partial coherence. 
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FIGURES 

FIG. 1. Tunneling in partial coherence through a series of barriers. The square represents a 

tunnel barrier, and the triangle stands for a phase-breaking scatterer which couples the system to 
a dephasing reservoir via two channels. 

FIG. 2. Schematic diagram for the scattering on the rath dephasing scatterer. The outgoing 
amplitudes are connected with the incoming amplitudes via = S^a^, see text for more details. 

FIG. 3. Additional resistance caused by momentum randomization on the dephasing scatterer. 
Here the special case of two identical resistors connected in partial coherence is demonstrated 
to show the interplay between dephasing and tunneling. 5 characterizes the tunneling strength 
through the individual resistor, i.e. Ti = T2 = 5. 

FIG. 4. Partial-coherence tunneling through N + 1 identical barriers (i.e. Tj = 6, 
j = 1,2,---,N + 1). The relative tunneling probabilities Teff (e)/reff(e = 0) and feg (e)/feff (e = 0) 
are respectively shown by the solid and dashed curves. The weak (a) and strong (b) tunneling 
limits are plotted to highlight the similarity and difference of the two dephasing models in different 
regimes. 
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